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ABSTRACT 



Context. The secondary-to-primary B/C ratio is widely used to study Galactic cosmic-ray propagation processes. The 
^H/*He and '^He/'^He ratios probe a different Z/A regime, therefore testing the 'universality' of propagation. 
Aims. We revisit the constraints on diffusion- model parameters set by the quartet (^H,'^H,'^He,^He), using the most 
recent data as well as updated formulae for the inelastic and production cross-sections. 

Methods. The analysis relies on the USINE propagation package and a Markov Chain Monte Carlo technique to estimate 
the probability density functions of the parameters. Simulated data are also used to validate analysis strategies. 
Results. The fragmentation of CNO cosmic rays (resp. NeMgSiFe) on the ISM during their propagation contributes to 
20% (resp. 20%) of the and 15% (resp. 10%) of the ^He flux at high energy. The C to Fe elements are also responsible 
for up to 10% of the ''He flux measured at 1 GeV/n. The analysis of '^He/'^He (and to a less extent ^H/*He) data shows 
that the transport parameters are consistent with those from the B/C analysis: the diffusion model with S ~ 0.7 
(diffusion slope), Vc ~ 20 km s~' (galactic wind), Va ~ 40 km s~' (reacceleration) is favoured, but the combination 
5 ~ 0.2, Vc ^ 0, and Va ~ 80 km s~' is a close second. The confidence intervals on the parameters show that the 
constraints set by the quartet data are competitive with those brought by the B/C data. These constraints are tighter 
when adding the '^He (or ^H) flux measurements, and the tightest when further adding the He flux. For the latter, the 
analysis of simulated and real data show an increased sensitivity to biases. Using secondary-to-primary ratio along with 
a loose prior on the source parameters is recommended to get the most robust constraints on the transport parameters. 
Conclusions. Light nuclei should be systematically considered in the analysis of transport parameters. They bring 
independent constraints which are competitive with those obtained from the B/C analysis. 

Key words. Astroparticle physics - Methods: statistical - ISM: cosmic rays 



1. Introduction 

Secondary species in Galactic cosmic rays (GCRs) are pro- 
duced during the CR journey from the acceleration sites 
to the solar neighbourhood, by means of nuclear inter- 
actions of heavier primary species with the interstellar 
medium. Hence, they are tracers o f the CR transport in 
the Galaxy (e.g.. IStrong et al]l2007D . Studying secondary- 
to-primary ratios is useful as it factors out the 'un- 
known' source spectrum of the progenitor, leav ing ^H/^He, 
^He/^He, B/C, sub-Fe/Fe— a nd recently p/p (jPutze et alJ 
l2009t Idi Bernardo et al1l2010D — suitable quantities to con- 
strain the transport parameters for species Z < 30. 

Most secondary-to-primary ratios have A/Z ~ 2, and in 
that respect, ^He/^Hc is unique since it probes a different 
regime and allows to address the issue of the 'universality' 
of propagation histories. For instanc e, in an analysis in the 
leaky-box model (LBM) framework, IWebberl (|1997[ ) found 
that '^Hc/*Hc data imply a similar propagation history for 
the light and heavier species (which was disputed in earlier 
papers). Webber also argued that the situation with regard 
to the ^H/^He ratio is less clear, because the uncertainties 
on the measurements are large (mainly due to instrumen- 
tal and atmospheric corrections). H and He spectra are the 
most abundant species in the cosmic radiation, and thus ^H 



and •^He are the most abundant secondary species in GCRs. 
However, achieving a good mass resolution — especially at 
high energy — is experimentally challenging. This explains 
why the elemental B/C ratio received more focus both ex- 
perimentally and theoretically (thanks to its higher preci- 
sion data w.r.t. to the quartet data). 

From the modelling side, after the first thor- 
ough and pioneering stud i es performed in the 60 ' s-70's 
(iBadhwar fc Daniel 19631: [Ramatv fc Lingenfeltei 1969 : 
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Meveil 119721: iMitleij 119721: iRamadurai fc Bis'^ Il974 : 
Mewaldt et al 1 119761) . the interest for the quartet nuclei 
somewhat stalled. Several updated analyses of the prop- 
agation parameters from the quartet were published as 
new data became available (see Table lA.ll for references) . 
How ever, very few dedicated studies w ere carried out in the 
80's (lBeattvlll986l: IWebber et al.lll987D likewise in the 90's 
(IWebbed Il990al : ISeo fc PtuskinI 11994 IWebbed fT99l . and 
none in the OO's. This is certainly related to the very slow 
pace at which new data became available in this period. 
Curiously, the most recent published data have not really 
been properly interpr eted, i.e. for '^ H/'^He data, IMAX9 2 
(Ide Nolfo et al.l [2000I) and AM S-Ol (lAguilar et all [20TTI): 
an d for ^He/^He dat a, IMAX92 ('Mcnn et al."2000(), SMILI- 
II (lAhlen et al.l2000ll . AMS-Ol (Xi ong ct al, 200ii), BESS9 8 
(iMvers et al.ll2003D . CAPRICE98 (iMocchiutti et al.ll2003D . 
Furthermore, almost all analyses have been performed in 
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the successful but simplistic LBM, but in a few studied. 
At the same time, the analysis of the B/C ratio has been 
scrutinised in more details. For instance, to replace the old 
usage of matching the data by means of an inefhcient man - 
ual scan of the parameter space (e.g., iJones et "all 1200 It) , 
more systematic scans were carried out (on the B/C and 
sub-Fe/Fe ratio) to get best-fit valu es as wc h as uncertain- 
ties on the parameters ( Maurin et al. 200 l i: Lionetto et al.l 
120051 lEvoh et al.l 120081: Idi Bernardo et al.ll2010D . A recent 



improvement is the use of Markov Chain Monte Carlo 
(MCMC) techniques to directly access the probability- 
density funct ion (PDF) of the GCR t r anspo r t and source 
param eters (jPutze et al.l l200l l2010l 1201 it iTrotta et all 
[20ll . 

In this paper, we revisit the constraints set by the quar- 
tet nuclei and their consistency with the results of heavier 
nuclei. In the context of the forthcoming PAMELA and 
AMS-02 data on these ratios, we also discuss the strategy 
to adopt and intrinsic limitations of the transport parame- 
ters reconstruction. For that purpose, we take advantage of 
the data taken in the last decade as well as simulated data 
of any precision, and analyse them with an MCMC tech- 
nique implemented in the USINE propagation code. This 
extend s and complement s analy ses of the B/C and primary 
nuclei (|Putze et al.ll2010l 1201 ll i in a ID diffusion model. 

The paper is organised as follows. In Sect. we briefly 
recall the main ingredients of the ID diffusion model and 
the MCMC analysis. We also list the parameters which are 
constrained. The simulated data and their analysis are de- 
scribed in Sect. [3] The analysis of the real data is given 
in Sect.m We conclude in Sect. [5] Appendix 1X1 gathers the 
data sets and the updated cross-sections used in the quartet 
analysis. 



2. MCMC technique, propagation and parameters 

The MCMC technique and its use i n the USINE propaga- 
tio n code is det a iled inlPutze et al.l (|2009f ) and summarised 
in iPutze et al.l (|2010t ). The ful l details regarding the ID 
transport model can be found in lPutze et al.l (|2010|) . Below, 
we only provide a brief description. 



2.1. An MCMC technique for the PDF of the parameters 

The MCMC method, based on Bayesian statistics, is used 
to estimate the full distribution (conditional PDF) given 
some experimental data and some prior density for these 
parameters. Our chains are based on the Metropolis- 
Hastings algorithm, which ensures that the distribution of 
the chain asymptotically tends to the target PDF. 

The chain analysis refers to the selection of a subset 
of points from the chains (to get a reliable estimate of the 
PDF). The steps at the beginning of the chain are discarded 
(burn-in length) if they are too far of the region of interest. 
Sets of independent samples are obtained by thinning the 
chain (over the correlation length). The final results of the 



99411 used a ID diffus i on mo del with reac- 



Webber fc Rockstrohl (|1997l ') relied on 



^ ISeo fc Ptuskiiil { 
celeration whereas 
Monte Carlo calculation; both studies conclude similarly (con- 
sistency with the grammage required for heavier species to 
produce the light secondaries). A preliminary effort based 
on the GALPROP pro pagation code was also carried out in 
iMoskalenko et al.l (|2003l ). 



MCMC analysis are the joint and marginalised PDFs. They 
are obtained by counting the number of samples within the 
related region of the parameter space. 

2.2. ID Propagation model and parameters 

The Galaxy is modelled to be an infinite thin disc of half- 
thickness /i, which contains the gas and the sources of CRs. 
The diffusive halo region (where the gas density is assumed 
to be equal to 0) extends to +L and —L above and below 
the disc. A constant wind V{r) ~ sign(z) -Vc x e^, perpen- 
dicular to the Galactic plane, is assumed. In this framework, 
CRs diffuse in the disc and in the halo independently of 
their position. Such semi-analytical models are faster than 
fuU numerical codes (GALPROFEI and DRAGOINEI), which 
is an advantage for sampling techniques like MCMC ap- 
proaches. 

2.2.1. Transport equation 

The differential density of the nucleus j is a function 
of the total energy E and the position r in the Galaxy. 
Assuming a steady state, the transport equation can be 
written in a compact form as 



_d_ 

'dE 



'dFj 



(1) 



The operator C (we omit the superscript j) describes the 
diffusion K{r,E) and the convection V{r) in the Galaxy, 
but also the decay rate ri.ad(^^) = 1/(7To) if the nu- 
cleus is radioactive, and the destruction rate Ti,^ei{r, E) = 
S/SA/ "'iSM('')''^crinoi(£') for collisions with the interstellar 
matter (ISM), in the form 

£(r, E) = -V- (KV) + V-V + F.ad + Ti^ei- (2) 

The coefficients b and c in Eq. ([T|) are respectively first 
and second order gains/losses in energy, with 



. , „N / dE \ 



v.v 



ion. coul. 



Ek 



2m + Ek 
m + Eh 



c{r,E) 
In Eq. 



(l + /3^) 
E 



X Kr, 



(3) 



(4) 



are 



— P ^ .[Vpp. 

([3]), the i onisation and C oulomb energy losses 
taken from Mannheirn fc Schlick eiscr (1994,) and 



IStrong fc Moskalenkol ()1998[) . The divergence of the 
Galactic wind V gives rise to an energy loss term re- 
lated to the adiabatic expansion of cosmic rays. The last 
term is a first order contribution in energy from reacceler- 
ation. Equation Q corresponds to a diffusion in momen- 
tum space, leading to an energy gain. The associated dif- 
fusion coefficient Kpp (in momentum space) is taken from 
the mod el of minimal reacceleration by the interstellar t ur- 
bulence (|Osborne fc Ptuskinlfigsl ISeo fc Ptuskinlll994D . It 
is related to the spatial diffusion coefficient K by 



A'pp X K 



3 " (5 (4 -(52) (4 



(5) 



where Va is the Alfvenic speed in th e medium. 

We refer the reader to App. A of lPutze et"all ()2010() for 
the solution to Eq. (P) in the ID geometry. 



http : //galprop . stanf o rd.edu/j 



http : //www ■ desy . de/ ~maccione/DRAGOH/ 1 



2 



B. Coste et al.: Constraining Galactic cosmic-ray parameters with Z <2 nuclei 
Table 1. Models tested in the paper. 



Model 


Transport parameters 


Description 


II 


{Ko, S, Va} 


Diffusion + reacceleration 


III 


{Ko, S, V,,Va} 


Diff. + conv. + reac. 



Note 1. For the sake of consistency, the model identification follows 
that of lPutze et"all I1200a . i20l"ol . 120111') and lMaurin et al.l ||2010|') . 



2.2.2. Free parameters of the analysis 

The exact energy dependence of the source and transport 
parameters is unknown, but they are expected to be power 
laws of 7?. = pc/Ze (rigidity of the particle). 

The low-energy diffusion coefficient requires a /3 = v/c 
factor that takes into account the inevitable effect of parti- 
cle velocity on the diffusion rate. However, the recent anal- 
ysis of the turbulence dissipation effects on the transport 
coefficient has shown that this coefficie nt could i ncrease at 
low-energ y ([Ptus kin ct al. 2006; Shalc hi fc Biisch ing 201^. 
Following iMaurin et al.l ((20100, it is parametrised to be 

KiE) = ■ KoTl^ . (6) 

The default value used for this analysis is 7/t = 1- The two 
other transport parameters are Vc, the constant convectivc 
wind perpendicular to the disc, and Va, the Alfvenic speed 
regulating the reacceleration strength [see Eq. ([S])] . The two 
models considered in this paper are given in Table [1] 

The low-ener gy prirnary source spectru m from accelera- 
tion models fe.g.. lDrurvlll983l : lJoneslll994|) is also unknown. 
We parametrise it to be 



dQ 



dE, 



k/n 



(7) 



where q is the normalisation. The reference low-energy 
shape corresponds to 77s = —1 (to have dQ/dp oc 
i.e. a pure power- law). 

The halo size of the Galaxy L cannot be solely deter- 
mined from secondary-to-primary stable ratios and requires 
a radioactive species to lift the degeneracy between Kq and 
L. However, the r ange of allowed va lues is still very loosely 
constrained (e.g.. lPutze et al.ll2010D . As the transport and 
source parameters can always be rescalcd would a differ- 
ent choice of L ass umed (see the scaling relations given in 
iMaurin et al.|[20Tol where S is shown not to depend on L), 
we fix it to L = 4 kpc. This will also ease the comparison of 
the resu lts obtained in this pape r with those of our previous 
studies (jPutze et al.ll2010ll201lD . 



3. MCMC analysis on artificial data sets 

MCMC techniques make the scan of high-dimensional pa- 
rameter spaces possible, such that a simultaneously es- 
tima tion of transport and source parameters is possi- 
ble (jPutze et al.l l2009t ). However, transport parameters 
are shown to be strongly degenerated for the B/C ra- 
tio d ata in the range 0.1 — 100 GeV/nuc (jMaurin et al.l 
I2010H. and s ource a nd tr ansport parameters are correlated 
(|Putze et al] l2009l I201GD . For GCR data in general, the 
fact that primary fluxes and secondary fluxes are not mea- 



sured to the same accuracjQ can bias or prevent an ac- 
curate determination of these parameters: a simultaneous 
fit has been observed to be driven by the more accurately 
measured primary flux (jPutze et al.ll201lh . This, although 
statistically correct, might not maximise the information 
obtained on the transport parameters. Therefore, several 
strategies can be considered when dealing with GCR data: 

— a combined analysis of secondary-to-primary ratio and 
primary flux to constrain simultaneously the source and 
transport parameters; 

— a secondary-to-primary ratio analysis only, either fixing 
the source parameters (i.e., using a strong prior), or 
using a loose prior. 

— a primary flux analysis only, either fixing the transport 
parameters (i.e., using a strong prior), or using a loose 
prior. 

In the literature, the strong prior approach has almost al- 
ways been used to determine the transport or the source 
parameters. The issue wc wish to address is how sensi- 
tive the sought parameters are to various strategies. This 
is the motivation to introduce artificial data, i.e. an ideal 
case study, as opposed to the case of real data where sev- 
eral other complications can arise (systematics in the data 
and/or the use of the incorrect propagation model or solar 
modulation model/level). 

3.1. Sets of artificial data 

To be as realistic as possible, we choose models that roughly 
reproduce the actual data points (see Fig. [4|), but also 
match the typical energy coverage, number of data points, 
central value and spread (error bars) of the measurement^. 
To speed-up the calculation and for this section only, we as- 
sume that all ■^He comes from ^He (see Sect. l4.11 for all the 
relevant progenitors). No systematic errors were added al- 
though they may set a fundamental limitation in recovering 
the cosmic-ray parameters. In practice, the statistical errors 
for the artificial data sets correspond to the sigma of the 
standard Gaussian deviations used to randomise the data 
points around their model value: "^Hc/^Hc was generated 
with statistical errors of 10% while He fluxes were generated 
with 1% and 10% errors, to simulate the situation where 
primary fluxes are 'more accurately' or 'equally' measured 
(in terms of statistics) than the secondary-to-primary ratio. 

The parameters of the two models used to simulate the 
data are listed in the two italic lines in Table [2l denoted 
Model II and Model III. They correspond to extreme values 
of the diffusion slope S, but which still roughly fall in the 
range of values fou nd for instance from the B/C analysis 
(|Putze et al.l |201G[ ) : for Model II with reacceleration only 
{Vc = 0), S is generally found to fall between 0.1 and 0.3, 
whereas for Model III with convection and reac celeration, S 
is generally found to fall i nto the 0.6—0.8 range (jJones et al.l 
1200 11 : IMaurin et al.ll201fll) . 



Statistical uncertainties are smaller for primary fluxes (more 
abundant than secondary fluxes), but the latter a more prone 
to systematics than ratios (e.g. secondary-to-primary ratios used 
to fit transport parameters). 

^ The uncertainty on the H and He fluxes is a few percents 
(for the recent PAMELA data, lAdriani et al.|[201ll ') and several 
tens of percents for the "^He/'^He ratio. 
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Table 2. Simulated data analysis for several Models (input parameters in italic) with L = 4 kpc: each line corresponds 
to the MCMC-reconstructed values (most-probable value, and relative uncertainties corresponding to the 68% CI) based 
on a given data/parameters option (see Sect. 13. 2p . The last column gives the value of the best x^/d-o.f. configuration 
found (corresponding to the curves shown in Fig. [1]). 



Option: data/params 


?7T 

- 


Ko X 10^ 
(kpc^ Myr"') 


5 Vc 

(kms-^) 


(kms~^) 


a 
- 


- 


Xbost/d.o.f. 
- 








Model 11 












1 


10.0 


0.2 


70 




1 




1: 3He/4He+He^o''° 
2: 3He/''He+He'''° 
2': 2 + src=true 
3: 3He/*He 
3': 3 -f src=true 


[1] 
[1] 
[1] 
[1] 
[1] 


I Q 3+3.9% 

10 2+2-0% 

Q 7+3.1% 
^•'-3.1% 

II r + 13.% 
11.0_23 % 

;^+3.0% 


102+3-1% 

19+21% 
19fi+5-i% 


-5.6% 
73 0+2.7% 

CO n+4-4% 

b».U_2 g% 

30 4+68.% 
'^^■*-53.% 

"^■"-6.6% 


9 9Q+0.9% 
^•^^^-1.3% 

9 90+0.3% 
^■^»_0.3% 

[2.3] 

2 70+148% 

^•'W_85% 

[2.3] 


0.78lf,^: 
[1] 

1 r+326% 

1-J_106% 
[1] 


0.93 
1.08 
0.97 
0.87 
0.88 


4: ^He/^He + src=prior 


[1] 


g Q + 11.% 
^•'-'-^11.% 


n 9+15.% 

^•^-10.% 


70 q+5.5% 

''-'■'^-8.2% 


[1.8,2.5] 


h2,+2] 


1.08 



Model III 







i.5 


0.75 


0.7 


18 


41 


2.3 


1 




1: 


3He/*He+He"''° 




1 r+40.% 


n CI +33.% 


18 /l + 17-% 


54 -1+18.% 


9 20+0-9% 
^•^^-1.7% 


1 +7.0% 

l^--36.% 


1.43 


2: 


3He/''He+Hei''° 




62+^ '^^° 


0.7251°:?^: 


10 7+3.6% 

1»-'_4.1% 


07 1+9.2% 
•-" -1^-1.3% 


9 004+0.3% 


A Q8+4-l% 

u.yo_2 0% 


1.00 


2' 


2 -f src=true 


1 AS+0-''% 
-l-48_o.7% 


7+4.3% 
"•'-5.7% 


o.nlill 




40 7+2-7% 
4U./_2 2% 


[2.3] 


[1] 


1.01 


3: 


3He/*He 


1 /i7+8.8% 
J^-*'-12.% 


Q 02+89 % 




2o.5ii«:^: 


58.01?^:^: 


12+^08% 

"•-1^-92.% 


2 2+^0.% 

■ —36.% 


0.83 


3' 


3 + src=true 


1 34+29.% 


n /IS+139% 

0-4o_6o.% 


0.68^22;^° 


18 -1+6.1% 
l^°-l^-3.3% 


44.8l^t.l 


[2.3] 


[1] 


0.89 


4; 


^He/'^He -f src=prior 


1 00+13.% 
l"JO_16.% 


pi 07+132% 
U.O(_g2.% 




9ri q+13.% 
^'-'••-'-4.0% 


42.21--;: 


[1.8,2.5] 


[-2, +2] 


0.96 



Model III: analysis with Model II 







1.5 


0.75 


0.7 


18 


41 


2.3 


1 




1: 


3He/*He+Heio''° 


[1] 


1 Q D + 3.6% 

J^'3-»_5.1% 


r. 91 +4.8% 
^J-^l_4.8% 




19C+4.8% 
1^D_4.0% 


9 0+0.9% 

^•'J_1.7% 


24+^2-% 

"•^*-29.% 


4.0 


2: 


3He/4He+He^^° 


[1] 


11-41?:^^: 


of;3+3-S% 




or + 2.4% 
°°-3.5% 


9 4+0.4% 
^•*-0.4% 


1 ig+1-7% 
-^■-^^-3.4% 


18 


4: 


^He/^'He + src=prior 


[1] 


20.71-2^: 


0.089l3i-% 




1071^2'% 


[1.8,2.5] 


[-2, +2] 


2.1 



Model II: analysis with Model III 







1 


10.0 


0.2 





70 


2.3 




1 




1: 


3He/*He+Heio''° 


1 07+5.1% 
-i-'J' -4.4% 


^•'-'-45% 


29+2'*''° 


1 7 9+9.3% 

l^'-^_41% 


70 r + 7.2% 


9 91+1.8% 
^•^-^-0.5% 


0. 


=0+12% 

'"-29% 


0.97 


2: 


3He/*He+He^''° 


r, 07+9.2% 


7 4 + 19% 
'•^-27% 


ri 9r+20% 
"•^O_8.0% 


f^+62% 

8-0-50% 


7/^ 0+3.5% 
'""^-4.8% 


9 90+0.4% 
^■^'-'-0.4% 


o.e 


0+7.4% 
'°-10% 


1.1 


4: 


^He/^'He + src=prior 


r, 40 + 193% 
"•^'-'-107% 


4 g+49% 

^•^-73% 


0.241-^: 


1 4+12% 
^8-4_5g% 


«g g+12% 


[1.8,2.5] 


[- 


-2, +2] 


1.05 



Note 2. A value in square brackets corresponds to the fixed value of the parameter for the analysis. An interval in square brackets corresponds 
to the prior used for the analysis (the posterior PDF obtained is close to the prior). 



3.2. Strategies to analyse the data 

To test the impact on the reconstruction of the transport 
{rjT, Kq, S, Va, and Vc) and/or source (a and rjs) param- 
eters, we test the following combinations (data set | model 
parameters) for the analysis. 



Option 4: transport (source = weak prior). 

We find that the He data alone cannot constrain the 
transport paramet ers (not shown here), in agreement with 
iPutze et al.l ()2011[ ) results (strong degeneracy between a 
and (5, but also with Kq, Va, and rjx). 



^HeZ-^He + He data 

Option 1 ((Thc = 10%): transport + source; 

Option 2 ((Thc = 1%): transport + source; 

Option 2' (crHe = l%): transport (source = 'true' value); 

^He/^'He ratio only 

Option 3: transport + source; 

Option 3': transport (source = 'true' value); 



3.3. Analysis of the artificial data 

In a first step, wc used the MCMC technique to estimate 
the best-fit parameters. The ■^He/^He ratio and ^Hc flux are 
shown for model II (crosses) and the corresponding simu- 
lated data (plusses) in Fig.[T] When both the ^He/''He and 
''He data are included in the fit (options 1 and 2, red dotted 
and red solid lines), the initial flux (crosses) is perfectly re- 
covered for ''He, and very well recovered for ■^He/^He. When 
the fit is only based on •^He/^He (options 3 and 5, magenta 
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E^O"[GeV/n] 



10" 



> 

o 



10' 



10- E 



P-I 

X 



' 1 1 1 

r "^He (simulated) x 


Input model - 


+ 


Simulated data (a=\%^^ 


I .*' X Best-fit from: 




yC / Option 1 




~ / / / Option 2 




=/ -J^ / — Option 3 




- X / — Option 3' 




- ^ / — Option 4 





10 



[GeV/n] 



Fig. 1. Analysis of simulated interstellar (IS) data sets for 
the ^He/"'He ratio (top panel) and the "^He flux times -E^y^ 
(bottom panel) based on Model II (x symbols), the pa- 
rameters of which are given in Table [2]) ■ The best- fit recon- 
structed curves correspond to the different 'options' given 
in Sect. [321 





Simulated data 

— Input paraiii. 

using ■'He/*He and He data: 

Option 1 (Ojy^ = 10%) 

— Option 2 (a4jj^ = 1%) 



or using He/ He data only: 

— Option 3 (src = free) 

Option 3' (src = true) 

— Option 4 (src = weak prior) 




Fig. 2. Marginalised posterior PDF for the transport and 
source parameters on the artificial data for Model II (the 
values for the input model are shown as thick vertical grey 
lines in each panel). The colour code and style correspond 
to the five 'options' described in Sect. 13.21 (and used in 
Fig. p. 



dashed and blue solid lines) , the initial flux is obviously not 
recovered (unless the source parameters are set to the true 
value as in option 3'), but the '^He/^'He ratio is consistent 
with the data. Unsurprisingly, the associated Xbcst/'^-^-^- 
values (last column of Tabled are close to 1. 

The MCMC analysis allows us to go further as it pro- 
vides the PDF of the parameters, from which the most- 
probable value and confidence intervals (CIs) are obtained. 
The results are gathered in Table [5] and Fig. [2] The vari- 
ous panels of the latter represent the PDFs for transport 
and source parameters for each 'option' for Model II. (For 
concision the correlation plots are not shown.) From these 
plots, some arguments are in favour of a simultaneous use of 
the secondary-to-primary ratio and the primary flux (here, 
^He/^'He and He), but not all. 

3.3.1. Advantages from a simultaneous analysis (ratio + flux) 

A simultaneous analysis (^He/''He 4- He) gives more strin- 
gent constraints on the transport parameters than only 
analysing the secondary-to-primary ratio (compare the 
PDFs for the red curves and blue curves in Fig. [2] respec- 
tively, for Ko, S, and 14). This partly comes from the ob- 
served correlations between transport and source param- 



etersB (e.g., iPutze et all l2009t) . Option 2' and option 3' 
with fixed source parameters show that the better CIs on 
the transport parameters come from the information con- 
tained in primary fluxcfO. The same conclusions hold true 
for Model HI (Vc 0), although with larger relative uncer- 
tainties due to the two extra transport parameters of the 
model {rjT and 14). 

3.3.2. What if the wrong model is used? 

As an illustration, we analyse data simulated from model 
HI {Vc 0) with model II {Vc = 0) and vice versa 
(lower half of Table If we force = (while = 
18 km s~^). the diffusion slope goes to a low value i5 ~ 0.2 
^^truo _ 0.7), while the Alfvenic speed goes to a high 
value Va ^ 100 km s^^ (V;'™'= = 41 km s^^). The larger 



" More stringent constraints on the source parameters (from 
more precise data) leads to more stringent constraints on the 
transport values: option 1 (crHe = 10%) vs option 2 ((thc = 1%). 

^ Note that the lack of constraints on the source parameters 
for option 4 confirms that the secondary-to-primary ratio is only 
marg inally sensitive to the source parameters (e.g.. IPutze et al] 

Eonl). 



5 



B. Coste et al.: Constraining Galactic cosmic-ray parameters with Z <2 nuclei 



X^(,g^/d.o.f. value with respect to the one obtained fitting 
the correct model easily disfavours this model. The second 
test (simulated with II, analysed with III) indicates whether 
allowing for more freedom in the analysis (two additional 
free parameters rjT and Vc) affects the recovery of the pa- 
rameters. The values of (5*™'= = 0.2 and = 70 km s'^ 
are recovered, while the others are systematically offset but 
less than Scr away from their true value. In this simple 
example, adding extra parameters is not an issue as the 
X^gs(./d.o.f. still favours the minimal model. H owever, with 
real d ata (see Sect.l4.3laii d the B/C analysis of lPutze et al.l 
[2oToH m aurin et al.l 120101) . in such a situation, it is so far 
impossible to conclude whether the correct model is used, 
due to the possible issue of multimodality and biases from 
systematics (see below). 



3.3.3. Drawbacks from a simultaneous analysis 

For Model II but even more for Model III (which has more 
free transport parameters), a possible worry is the pres- 
ence of multimodal PDF distributions, which more often 
happens for the simultaneous analysis. An example of mul- 
timodality is the analysis with Model II (i.e. Va = 0) of 
data simulated with Model III, which corresponds to a lo- 
cal minimum of the true Model III parameters. This is of 
no consequence for the ideal case, but real data may suf- 
fer from systematics errors and/or the inappropriate solar 
modulation model may be chosen. In that case, the true 
minimum can be displaced, or turned into a local mini- 
mum (and vice versa) . Measurements over the last decades 
showed that primary fluxes are more prone to systemat- 
ics than secondary-to-primary ratios. Primary fluxes are 
also more sensitive to solar modulation than ratios. For 
these reasons, the use of secondary-to-primary data only 
(option 4) for the analysis, although less performant to get 
stringent limits on the transport parameters, is expected to 
be more reliable and robust. 



Obviously, our analysis does not cover the range of all 
systematics when dealing with real data. A more system- 
atic analysis — e.g. covering a wider family of propagation 
models, several solar modulation models, several sources 
of systematics in the data — goes far beyond the scope of 
this paper. Note that some of these effects are likely to be 
energy dependent, complicating even further the analysis. 
With the successful installation of the AMS-02 detector on 
the ISS and its expected high-precision data, these issues 
arc bound to gain importance. 



4. Constraints from the quartet data 

We now apply the MCMC technique to the analysis of real 
data. We emphasise that for the artificial ones, we assumed 
the ■^He to come solely from the *He fragmentation, in order 
to speed up the calculation. Based on our new compilation 
for the cross-section formulae (see App. |B]), we take into 
account the contributions from A > A CR parents, checking 
which parents are relevant ( tj4.1[) . Having determined the 
heaviest parent to consider in the calculation, we then move 
on to the result of the MCMC analysis ( tj4.2p . and those 
from our best analysis fi j4.3p . 



4.1. Fractional contributions 

At first order, the contribution to the ^H and "^Hc secondary 
production from Z > 4 nuclei is proportional to the source 
term [see Eq. ([1])]. For a secondary contribution, the 
source term is proportional to the primary flux of the par- 
ents (which have been measured by many experiments), and 
to the production cross-section. Normalised to the produc- 
tion from ^He, we have 



■7l, 



(8) 



3.4. Recommended strategy to analyse real data 

The most robust approach to determine the transport pa- 
rameters (and their CIs) is to analyse the sccondary-to- 
primary ratio using a loose but physically-motivated prior 
on the source parameters (option 4). This has the advan- 
tage of taking into account the correlations between the 
source and transport parameters. The simultaneous anal- 
ysis is mandatory to obtain the source parameters (op- 
tion 1 or 2). It also brings more information on the trans- 
port parameters, but the primary fluxes can bias their de- 
termination if it suffers from systematics. We recommend 
such an analysis to be performed in addition to the direct 
secondary-to-primary ratio analysis, in order to get the fol- 
lowing diagnosis: if the range of values for the transport 
parameters from both analyses are 

— inconsistent, it indicates that the values and CIs ob- 
tained for the sources parameters are biased or unreli- 
able; 

— consistent, the selected propagation model may be the 
correct one, and the source parameters are then the 
most probable ones for this model. However, the CIs on 
the transport parameters are very likely to be underes- 
timated if the error bars on the ratio are much larger 
than the ones on the primary fluxes. 



where P is the CR projectile, S is the secondary fragment 
considered, and 7I [see Eq. (|B.3|) ] is the production cross- 
section relative to the production from ^Hc. The fractional 
contribution /^^^ for each parent is defined to be 



Rcl 



P-i-S 



Ep'=Hc - Ni ^"^^ 



P'-s-S • 



(9) 



As seen from Tabic [3j the most important contributions 
from primary species heavier than He {Z > 2) are C and O, 
followed by Mg and Si and finally Fe. The total contribution 
of these species amounts to ~ 35% for ^H and ~ 11% for 
•^He, but mixed species (such as N) or less abundant species 
also contribute to ^ 5%. 

A proper calculation of these fractional contributions in- 
volves the full solution of the propagation equation, taking 
into account energy gains and losses, total inelastic reac- 
tions, and convection. Based on the propagation parame- 
ters found to fit best the current data (see next section), 
we show in Fig. |3] the fractional contribution oi A > A nu- 
clei as a function of energy for the full calculation. It con- 
firms the previous figures, but with a residual energy de- 
pendence (itself depending on the species) hitting a plateau 
above ^ 100 GeV/n. The difference can be mostly at- 
tributed to a preferential destruction of heavier nuclei at 
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Table 3. Estimated fractional contribution of projectile 
^ > 4 to the and ^He fluxes. The columns are respec- 
tively the name of the element, atomic number, ratio of 
measured fluxes, production cross-section ratios for and 
•^He, and the estimated fractional contribution to and 
3He. 



p 


Ap 




7p" 


1> 


07 
70 


(■P^-^He 

0/ 
70 


He 


4 


1.0 






60.1 


86.0 


C 


12 


3.310"^ 


5.5 


2.3 


7.5 


3.6 


N 


14 


7.410"^ 


6.6 


2.4 


2.0 


0.8 





16 


3.410"^ 


7.8 


2.6 


10.7 


4.1 


F 


19 


5.210"* 


9.6 


2.8 


0.2 


0.1 


Ne 


22 


5.1 10"^ 


11.4 


3.1 


2.4 


0.8 


Na 


23 


8.810"* 


12.1 


3.2 


0.4 


0.1 


Mg 


24 


6.710"^ 


12.7 


3.3 


3.4 


1.0 


Al 


26 


1.1 10"^ 


14.0 


3.5 


0.6 


0.2 


Si 


28 


5.510"^ 


15.3 


3.7 


3.4 


1.0 


P 


31 


1.810"* 


17.3 


4.0 


0.1 


< 0.1 


s 


32 


1.010"^ 


17.9 


4.1 


0.7 


0.2 


CI 


35 


1.810"* 


20.0 


4.5 


0.1 


< 0.1 


Ar 


36 


oil n — 'I 
3.1 10 


on c; 


4.D 


U.3 


n 1 
0.1 


K 


39 


2.1 10"* 


22.7 


5.0 


0.2 


< 0.1 


Ca 


40 


6.1 10"* 


23.4 


5.1 


0.6 


0.1 


Sc 


45 


1.010"* 


27.0 


5.8 


0.1 


< 0.1 


Ti 


48 


3.410"* 


29.2 


6.2 


0.4 


0.1 


V 


51 


1.810"* 


31.4 


6.6 


0.2 


0.1 


Cr 


52 


3.610"* 


32.1 


6.8 


0.5 


0.1 


Mn 


55 


3.010"* 


34.3 


7.2 


0.4 


0.1 


Fe 


56 


3.710"^ 


35.1 


7.4 


5.3 


1.3 


Ni 


58 


2.410"* 


36.6 


7.7 


0.4 


0.1 



Note 3. The ratio of measured fluxes is calculated at ~ 10 
GeV/n using PAMELA for He (lAdriani et al.ll201lD and HEAO-3 
llEngelmann et al.|[l990h for 6 < Z < 30. 



low energjlfl. Note that for production, the coalescence 
of two protons (long-dashed pink curve) contributes up to 
40% of the total at ~ 1 GeV/n energy (peak of the cross- 
section, see Fig. IB.3| ). Depending on the precision reached 
for the d ata, it is important to include the CNO contribu- 
tion (e.g. Ramatv fc Lingenfelteil Il969t lJung et aD Il973bt 
lBeattvlll986[ ). but also all contributions up to Ni. 

Finally, the fragmentation of CNO can also affect the 
and *He primary fluxes. The peak of contribution occurs at 
GeV/n as secondary fluxes drop faster than primary fluxes 
with energy. Fig.[2]shows this contribution to be < 10% for 
^He. With the high precision measurement from PAMELA 
and the even better measurements awaited from AMS-02, 
this will need to be further looked into in the future. 



Indeed, the primary-to-primary ratios are not constant. The 
heavier the nucleus, the larger its destruction cross-section, the 
more the propagated flux is affected/decreased at low energy, 
the longer it takes to reach a plateau of maximal contribution at 
high energy. The observed trend is co nsistent with the p rimary- 
to-primary ratios shown in Fig. 14 of iPutze et al.l (|201ll ). 



S 0.8 
o 



.2 0.6 



0.4 



0.2 - 





' — ' — 


'H,'He...O...Si...Zn \ 




- 'H,'He...O...Si 


- 


- 'H.'He .-5'\\ 


^ CR' + H/He^-H 1 

i— parents 


-: 







^10-' 



10 



10^ 

Ey„ [GeV/n] 



iP- 



i!0.95- 

"rt Q CR parents: 

■a : ■'He...O...Si...Zn 



1 0.85; 



0.8- 



0.75- 



■'He...O...Si 



^ CR' + H/He ^ -''He 

-parenls 



0.7. 



10-' 



10 



10^ 

E^„ [GeV/n] 



o 

■■=0, 



96 



0.94 





_ CR cohttibutions: 


^ cr' + H/He ^ ■'He - 


Prim + Li...O...Si...Zn 


i-parents /' 


PrimV Li...O...Si 


,/■' + 


_ Prim +\i...O 

4He (primafy) 


Primary 


1 1 





10-' 



10 



10" 

Ew„ [GeV/n] 



Fig. 3. Fractional contributions to the propagated 
fluxes (top panel), '^He fluxes (middle panel) and ''He fluxes 
(bottom panel) as a function of i?fc/„ from A> A CR par- 
ents. For ^He, the primary contribution is also considered. 



4.2. MCMC analysis: test of several data combinations 

Given the accuracy of current data (see Fig. H]), we must 
take into account the contribution from all parent nuclei at 
least up to ^ °Si. In the rest of the analysis, we use PAMELA 
data for He (jAdriani et al.l[20Tl[) , as they overcome all oth- 
ers in the ~ GeV— TeV range in terms of precision. Before 
giving our final results, and to complement Sect. 13.21 we 
discuss the appropriate choice of data to consider here, in 
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Table 4. MCMC analysis of Model III [Vc 7^ 0) with i = 4 kpc: most-probable values and relative uncertainties 
(corresponding to the 68% CI) for the analysis of various combinations of '^He/'^IIe, ^He, and PAMELA He data (-1- 3 
combinations involving the isotope). The last column gives the x^/d.o.f. for the best-fit model found (which usually 
differs from the most-probable one). 



Data 


Ko X 10" 


5 




Va 


a 


ns 


Xbcst/d.o.f. 




(kpc^ Myr^^) 




(kms-^) 


(kms^^) 









U.0U_j2% 

1 10+2''^'' 



^He/''He + He 
^He/^He + He* 



'He/'^He 

^He/^'He + ^He 

3He/*He ^He^ 

'He/^'He + ^He + He 

'He/^'He + ^He^ + He 



^H/^He 

^HZ-'He + ^H + He 
^H/^'He + ^H + He* 



1.20+^ 



r| + 58% 
^-8.3% 
r, 09+35% 
^•°^-50% 

"•'"-9.2% 

"•«^-8.6% 

1 r + 13% 
^^-33% 

0.53+jg5,° 

c p+55% 
"•°-45% 



0.79 
0.54 
0.67 

0.56 
0.39 



+6.4% 
5.2% 
+5.6% 
7.4% 
+6.0% 

15% 

+ 11.% 
8.9% 
+26.% 
5.1% 



n 70+3.8% 

"• '^-2.5% 
ri '7'7+5.2% 
U.( I _2S% 



o.os; 

0.61 
0.039 



100% 

21% 
8.2% 

6.6% 
+ 125% 
125% 



1 7 0+1.7% 
ll.O_2 5% 

-,+3.6% 



28.0] 
27.4 



3.6% 
+3.3% 
-4.7% 



24.2 



+4.5% 



-2.9% 
24 4+6-1% 

^+1.4% 
-2.8% 
-1.6% 
-3.1% 



19.6: 

19.2"* 



3.0 



+500% 
-67% 



14.6 
26.1 



+3.4% 
-2.7% 

+8% 
-27% 



39.5 
54.0 
41.0 



+6.1% 



+5.6% 
-9.3% 

+34% 

19% 



68.1 
85.0 



+6.8% 
16.% 
+7.1% 
22.% 
1 +4.0% 
"--4.8% 
5+5.0% 
"-3.7% 



33.7 
26.3 
24.1 



+52% 
28% 
+8.0% 
11% 
+39% 
46% 



1.8% 
2.7% 
+0.4% 
0.8% 



2.26 
2.49 
[1.8,2.5] 



9 90+2.7% 

^•^■J_3.6% 

9 -1 y+2.3% 
^•-L'_4.6% 

9 99+0.9% 
^•^^-1.8% 
9 91 +1.4% 
^•^1--1.4% 



2.38 
2.47 
2.741} 



+4.6% 
3.8% 
+1.2% 
2.4% 
11% 
1% 



0.07 



+59% 
74% 



-1.97 



+9.6% 
-1.0% 



[-2, +2] 



-0 

-0.82 
0.02 
0.04 



33% 

52% 
+ 18% 
34% 
171% 
159% 
+71% 
63% 



-1 1^^+17°''° 



0.47 
0.38 



+ 13% 
21% 
+47% 
39% 



4.8 

2.1 
2.9 

1.8 
1.9 
3.3 
3.4 

6.8 
6.4 
4.9 



^ Including AMS-01 data from the recently published analysis of lAguilar et al.l (|201ll ') 
* Excluding PAMELA He data points below 5 GeV/n and above 183 GeV/n. 



Note 4- An interval in square brackets corresponds to the prior used for the analysis (the posterior PDF obtained is close to the prior). 



order to get the best balance between robustness and reli- 
ability for the and ■^Hc-rclatcd analyses. 



4.2.1. Simulated vs. real data 



We start by comparing the results obtained with the simu- 
lated and the actual data set. To avoid lengthy comparisons 
of numbers, we limit ourselves to Model III (where we also 
fix r]T to its default value, i.e. 1). The obvious difference 
with the simulated data is that we no longer have access 
to the true source parameters (automatically excluding op- 
tions 2' and 4 discussed in Sect. I3.2)) . For the simultaneous 
analysis using He PAMELA data — the precision of which is 
~ 1% — , we recover similar values and CIs for the parame- 
ters (compare option 2 in Tableland the first three lines 
of Table |4]) . The second row of Table |4] is based on a subset 
of He data: high energy data points are discard ed because 
they show departure from a single power-law (jAhn et al.l 
l2010HAdriani et al.ll20Tl[) . whereas low -energy data points 
are discarded because of their sensitivity to solar modu- 
lation, which is presumably too crudely described by the 
Force-Field approximation used here. The x^j,g^/d.o.f. value 
(first row) shows that the model has difficulty to perfectly 
match the high precision PAMELA He data over the whole 
energy range. The analysis of '^He/^He ratio using a prior 
on the source parameters (option 4 in Tablc[2]and third line 
of Table H]) gives larger CIs for the transport parameters. 
The consistency between the results of the latter analysis 
(third line) and that based on the partial He data (second 
line) , and their discrepancy with the results of the analysis 
based on the full He data set (first line) confirms our sus- 
picion that high precision measurements for primary fluxes 
can bias the transport parameters determination. 



4.2.2. Adding the secondary '^Hc flux in the analysis 

Replacing He by ■^Hc in the simultaneous analysis (4"^ and 
5*^ line) further affects the determination of the transport 
parameters. This is not surprising since "^Hc data arc not 
all consistent with one another (see Fig. The bias is 
stronger when taking into accoun t the recently published 
AMS-01 data (jAguilar et al.ll201lt) . If both ^Hc and He are 
taken into accouni0, the much better accuracy of the He 
PAMELA data with respect to the '^Hc data amounts to a 
smaller weight of the latter in the analysis. 

4.2.3. ^H/^Hc vs ^Hc/^He 

We repeat partially the analysis for ^H in the last 3 rows 
of Table [21 The data arc so inconsistent with one another 
for ^H/''He (see Fig. 2]) that we are forced to use at least 
the ^H flux (whose data points are also markedly incon- 
sistent with one another). Even so, the results are not re- 
liable. PAMELA and AMS-02 have the capability to im- 
prove greatly the situation, but in the meantime, we are 
forced to include He as well in the analysis (next-to-last 
row in the table) . The transport parameter values from the 
^H/"'Hc+^H+Hc analysis are grossly consistent with those 
from the ■^He/*He-|-'^He-|-He analysis, but are likely to suffer 
from similar biases (see the previous paragraph). Reducing 
the energy range of He data is not even possible for the ^H 
analysis (last line in the table), as the results obtained are 
not reliable. 



^ The simultaneous analysis of '^He/*He + '^He -I- He has not 
been tested in the simulated data since it would have amounted 
to a double-counting of the '^He data (appearing in the three 
quantities). However, real data involve different experiments for 
the various quantities (PAMELA for He and other experiments 
for '^He), and independent measurements are used. 
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Table 5. Most-probable values and CIs for Models II and III for our 'best' analysis (Sect. 14. 3p of the quartet and B/C 
data (L = 4 kpc). 



Data 


Ko X 10^ 
(kpc^ Myr"^) 


S K 

(kms-i) 


Va 

(km s"^) 


a 


Vs 


Xbost/d.o.f. 


3jje/'*He 

^He/^'He + ^He + He* 
^H/^'He -f -f He 
B/C [Putze et al., 2010] 
B/C [this paper] 
B/C -f C (all) [this paper] 
B/C + C (HEAO) [this paper] 


1 c n+O.S 

14.8l°-j 
n 7+0-3 

^- ' -0.2 


Mode/ J/ 

n 2Q+0 02 
u.^y_o.o3 

o.3ii°:°i 

0.08lg°^3 

Q 2O4 + O.OO6 

""JJ-O.OI 

U.,314_Q QQg 

n qc;q+0.004 

U.OOO_Q gQ4 


116+" 
74+* 
44t| 
73t^ 
80i^ 
571? 
78tl 


fl 8 9 "il 

[i.O, Z.OJ 
9 qR+0 03 

2 66+° °^ 

a + S = 2.65 

[1.8,2.5] 
9 040+0.005 

9 ot;n+0 0i5 

Z.ZOU_g ggg 


7n+o-05 

-1 

[-2, +2] 

U.»D_g g4 

1 /IO+0.08 


0.0 

4.6 
5.6 
4.7 
1.5 
13.9 
2.8 



Model III 



^'He/^He 


n 5+''--* 

U-O-0.3 


n 67+0 04 
"•O'-O.IO 


27.4t?:° 


411^^ 


[1.8,2.5] 


h2, +2] 


2.9 


^He/^'He + ^He + He* 




0.64t°°^ 


23.51?:;; 


541^3 


9 07+0.03 

^"J' -0.03 


03+° °* 


1.6 


^H/^'He + He 


q n+0.8 
•J-^-0.1 


0.50l°°« 


27.ll?-| 


72t?g 


9 41+0.12 


sq+° °'' 


5.5 


B/C [Putze et al., 2010] 


46+° °® 

U-*O-0.06 




18.9l°i 


38l^ 


a -f 5 = 2.65 


-1 


1.5 


B/C [this paper] 


46+°-^® 

'-'■^"-0.10 


0.82«°« 


18.31°:^ 


40l^ 


[1.8,2.5] 


[-2, +2] 


0.9 


B/C + C (all) [this paper] 


'=i7+°-°5 
'-'•'J ' -0.03 


0.801°:°? 


17.41°:^ 


3611 


9 9Rn+° °°^ 

Z.ZUU_g ggg 


94+° °3 


5.2 


B/C -f C (HEAO) [this paper] 


r, qq+0.06 
'-'■'J '5-0. 10 


0.93l°°? 


i8.2i°:i 


351^ 


9 019+0.019 
Z.OiZ_g ggg 


1 q+0.1 


2.0 



Excluding PAMELA He point below 5 GeV/n and above 183 GeV/n. 



Note 5. An interval in square brackets corresponds to the prior used for the analysis (the posterior PDF obtained is close to the prior). 

Note 6. The B /C results are based on IMP7-8 jGarcia-Munoz et al.'' 1987l'l . Vo yager lfc2 lILukasiak et al.|[l999lh ACE-CRIS lIGeorge et al.l 
I200gfl . HEAO-3 llEngelmann et al.lll990fl . Spacelab {K/Iucllcr ct al. 1991), AMS-01 llAguilar et al.ll2011lh and CREAM llAhn et al.ll2008l'l . shown 
to be the most compatible data for a B/C analysis llPutze et al.ll2009l ). 



4.3. MCMC analysis: 'best' results 

Taking into account specificities of the actual data (previ- 
ous section), our 'best' analysis is based on the most rele- 
vant combinations of data for and ■^He: 



— the ^He/^He analysis (with a prior for the source 
parameters) gives robust and conservative results 
for the transport parameters. The result from the 
'^He/^He-l-'^He-l-He-'- analysis is more sensitive to biases, 
but using an energy sub-range for He data is expected 
to limit them. 

- Due to the paucity of data, the ^H/'^He-h^H-f He (full 
energy-range for He) analysis is the only reliable option, 
although it probably suffers from biases. 



The corresponding most-probable values and CIs are gath- 
ered in Table [5l and the corresponding envelopes for 
^H/^He, ^H, ^He/^He, and ^He/^He are given in Fig. H 
We also re-analyse the B/C ratio according to our 'best- 
analysis' scheme (B/C alone with a prior for the source pa- 
rameters or B/C + C). The re sults are rep o rted i n Table [5l 
where the results obtained in iPutze et al.l (j2010r ) for fixed 
source parameters are also reproduced: we note that the 
new strategy gives results in better agreement with those 
of the quartet analysis (e.g., the transport parameters 6 
and Kq are shifted by more than 30% for Model II), fur- 
ther demonstrating its usefulness. 



4.3.1. Universality of the transport parameters 

If we focus on the transport parameters, we note that com- 
binations involving the ^H/^He, '^He/^lie, or B/C ratio 
give broadly consistent transport parameter values, be it 
for Model II or Model 110. Regardless of the actual prop- 
agation model, we conclude that these results hint at the 
universality of CR transport for all species. Another impor- 
tant result is that the constraints set by the quartet data on 
the transport parameters are competitive with those set by 
the B/C ratio, so that the quartet data should be a prime 
target for AMS-02. 

4.3.2. Model II {S - 0.3) or Model III [S - 0.7)? 

According to Sect. 13.21 comparing the results of the 
secondary-to-primary ratio analysis with those of the com- 
bined analysis (ratio -f- primary flux) gives an indication 
of their robustness. Table [5] show that the results for the 
diffusion slope S is very robust, regardless of the model 
considered. A more detailed comparison shows that for 
■^He-related constraints, the transport parameter values for 
Model II arc inconsistent with one another at the Scr level, 
whereas the 68% CIs overlap with one another (but for Vc) 



^° The most significant difference is for the ^H/^He-l-^H-l-He 
analysis, which is inconsistent in both models and clearly unre- 
hable for Model II (5 ~ 0). For model HI, B/C and ^He/*He- 
related constraints are roughly in the same region but are lo- 
cated at several a from each others (they are consistent with 
one another for Model II). 
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I ' — ^ 

IMP3+4+5 (650 MV) 
IMP7+RonneerlO (540 MVT 
IMPS (400 MV) 
ISEE3 (740 MV) 
MASS89(12()0MV) 
IMAX92 (550 MV) 
BESS93 

AMS-01 (600 MV) 



Bkilloon73 (650 MV) 
Balloon77+Voyager (^ 
B!illoon84 (650 MV) 
Blllloon90(12()0MV: 

__i 



: -'He/*He (TOA) 



2 0.; 

J 



Balloon72 (400 MV) 
Balloon73 (500 MV) 
Balloon77+Voyager (400 MV) 
Balloon89(1400MV) 
IMP7+Pioneerl0 (540 MV) 
ISEE3 (740 MV) 
ISEE3(HIST) (500 MV) 



E^^ [GcV/nl 



Voyager87 (360 MV) 
IMAX92 (750 MV) 
SMILI-I(I200MV) 
SMILI-II (1200 MV) 
AMS-01 (600MV) 
BESS98 (700 MV) 
CAPRICE98 (700 MV) 





Model 11 wilh prior (^l = 700 MVl 
Model 111 wilh prior it) = 700 MV 
Model 11 ill = 700 MV) 
Model 111 (((. = 700 MV) 



Fig. 4. Left panels: demodulated interstellar (top) and ■^He (bottom) envelopes at 95% CIs times E'^J^_^. Right panels: 

top-of-atmosphere secondary-to-primary ratio ^H/'^He (top) and '^He/^He (bottom) ratios. The full envelopes correspond 
to the 'best' simultaneous analysis (secondary-to-primary ratio -I- primary flux) for model II (blue) and III (red). The 
hatched envelopes correspond to the '^He/'^Hc analysis (prior on source parameters). See Table ETTl for references and 
the corresponding demodulation (for IS) and modulation (TO A) level (j). 



for Model III, hence slightly favouring the latter [5 ^ 0.7). 
The comparison of the Xbost/'^-'^-^- values also tends to 
favour model III. Hence, although the value 5 ~ 0.7 seems 
favoured, we cannot exclude yet pure reaccel eration model 
Vr. = 0) with 5 ^ 0.3. Moreover, as shown in lMaurin et al.l 
2010f ). many ingredients of the propagation models can 
lead to a systematic scatter of the transport parameters 
larger than the width of their CIs. Data at higher energy 
for any secondary-to-primary ratio are mandatory to con- 
clude on this issue. 



4.3.3. Source spectrum 

The present analy sis is more general than that used in 
iPutze et al.l ()2011f ). where the transport parameters were 
fixed. Although it is not the main focus of this paper, we 
remark that the values of the source slope a fr om the B/C 
+ C a nalysis are consistent with those found in lPutze et al.l 
(|201lh . strengthening the case of a universal source slope 
a at the ~ 5% l evel. For the quartet values, ano is broadly 
consistent with iPutze et ahl 's analysis (based on AMS-01, 
BESS98 and BESS-TeV data for He). However, the results 
for the source parameters depend on the choice of data sets 
and energy-range considered. This indicates that for < 1% 
accuracy data, either the model for the source is inappro- 
priate, or the solar modulation model is faulty, or some 
systematics exist in the measurements. The AMS-02 data 
will help to clarify this question. 



5. Conclusion 

We have revisited the constraints set on the transport (and 
also the source) parameters by the quartet data, i.e. ^H, ^H, 
■^He, and *Hc fluxes, but also the sccondary-to-primary ra- 
tios ^H/'^He and ^ He/^He. This extends and complements 
a series of studies (jPutze et al.l l2009l . l2010l . 120111) carried 
out with the USINE propagation code and an MCMC al- 
gorithm. The three main ingredients on which the analysis 
rests are: 

— A minute compilation of the existing quartet data and 
survey of the literature, showing that the most re- 
cent/precise data (AMS-01, BESS93->98, CAPRICE98, 
IMAX92, and SMILI-II) have not been considered be- 
fore this analysis. 

— We have done a systematic survey of the literature for 
the cross-sections involved in the production/survival 
of the quartet nuclei. This has lead us to propose new 
empirical production cross-sections of ^H, ■^H, and '^He, 
valid above a few tens of MeV/n, for any projectile on 
p and He (we also updated inelastic cross-sections). 

— We have made an extensive use of artificial data sets 
to assess the reliability of the derived CIs of the GCR 
transport and source parameters for various combina- 
tions of data/parameters analyses. 

In broad agreeme nt wi th previous studies, (e.g. 

iRamatv fc Lingenfelteil Il969l : iBeattvl Il986f ). we find that 
the fragmentation of CNO contributes significantly to the 
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flux 30%) above a few GeV/n energies ("^He fragmen- 
tation is the dominant channel for the and '^He fluxes). 
Nevertheless, we provide a much finer picture, showing in 
particular that heavy nuclei (8 < Z < 30) contribute up to 
10% for ^He (20% for ^H) at high energy. We also provide an 
estimate of the secondary fraction to the ^He flux. By def- 
inition, the secondary contribution has a steeper spectrum 
than the primary one and therefore becomes quickly negli- 
gible at high energy. This secondary contribution peaks at 
a few GeV/n, and it amounts to ~ 10% of the total flux 
(~ 7% up to O fragmentation, ^ 2% from elements heavier 
than O), which is already a sizeable amount g iven the ~ 1% 
preci sion reached by the PAMELA data (jAdriani et al.l 
1201 If ). For ^H, the knowledge of the multiplicity of neu- 
tron and proton produced by the interaction of all elements 
on the ISM is required to calculate precisely its secondary 
content. 

Simulated data have allowed us to check several critical 
behaviours. Firstly, the He flux is obviously useful (and re- 
quired) to constrain the source parameters, but it has also 
been found to bring significant information on the transport 
parameters: fitting a secondary-to-primary ratio plus a pri- 
mary flux brings more constraints than just fitting the ratio 
(even when source parameters are fixed). Secondly, we have 
checked that a model with more free parameters (than the 
ones used to simulate the data) is able to recover the correct 
values. However, our analysis has also strongly hinted at the 
fact that adding the primary fiux He biases the determina- 
tion of the transport parameters if systematics (which are 
usually more important in primary fluxes than in ratios) are 
present, and/or if the wrong model is used. For this reason, 
when dealing with measurements, we recommend to always 
compare the result from the secondary-to-primary ratio + 
primary flux analysis to that of the secondary-to-primary 
ratio using a loose but physically motivated prior on the 
source parameters. 

The analysis of real data has shown that quartet data 
slightly favours a model with large 5 ~ 0.7 (with Vc ~ 
20 km s^^ and Va ^ 40 km s~^), but that a model with 
smaU S 0.2 (with K and K 80 km s'^) can- 
not be completely ruled out. Better quality data, and espe- 
cially data at higher energy are required to go further. The 
conclusions are similar and the range of transport parame- 
ters foun d are consistent with those obtained from the B/C 
analy sis (| Jones et al.ll200ll:lPutze et al.l[20lHMaurin etaD 
l2010fP^ . This strongly hints at the the universality of the 
GCR transport for any all nuclei. Furthermore, we have 
shown that the analysis of the light isotopes (and the al- 
ready very good precision on He) is as constraining as the 
B/C analysis (similar range of CIs). 

The several difficulties which have been pointed out in 
this analysis could be alleviated by virtue of using bet- 
ter data. However, it is more likely that the interpretation 
of future high-precision data will require the development 
of refined models for the source spectra and/or transport 
and/or solar modulation. For instance, the Force-Field ap- 
proximation for solar modulation is already too crude to 
minutely match the PAMELA He data. The forthcoming 



Note that in this paper, we did not attempt to combine 
the results of different secondary-to-primary ratios (^H/'*He, 
^He/''He, B/C, sub-Fe/Fe, p/p). This is left for a future study, 
for which a Bayesian evidence could be used to better address 
(in a Bayesian framework) the crucial issue of model selection. 



AMS-02 data at an even better accuracy will definitively 
pose interesting new challenges. 

Acknowledgements. D. M. would like to thank A. V. Blinov for kindly 
providing copies of several of his articles. A. P. is grateful for financial 
support from the Swedish Research Council (VR) through the Oskar 
Klein Centre. 



Appendix A: Cosmic-ray data 

Deuterons and '^He fluxes are very sensitive to the mod- 
ulation level, whereas ratios are less affected. The exact 
value for the solar modulation level </> is uncertain. For in- 
stance, the values given in the s eminal papers can differ 
greatly from those estimated by ICasadei fc Bindil ()2004|) 
(in order to match the electron and positron fluxes of var- 
ious experiments, see their Table 1), or from those recon - 
structcd from the Neutron monitors (jUsoskin et al.|[2"002[ ). 
This difference may be attributed to the fact that the lat- 
ter analysis correctly solves the Fokker-Plank equation of 
GCR transport in the Heliosphere, whereas most papers 
rely on the widely used force-field approximation that is 
know n to fail for strong modulation level > 1000 MV 
(e.g.. lUsoskin et al]|2002[ ). In this analysis, we do not at- 
tempt to go beyond this force-field approximation, as speed 
is of essence for our MCMC analysis. We rely mostly on 
the force-field effective modulation parameter (/) necessary 
to reproduce the data (as quoted in the seminal paper), but 
these values are slightly adjusted in order to give overlap- 
ping fluxes when all the data are demodulated and plotted 
together. Given the uncertainty on the data, the large un- 
certainty on 0, and the fact that most-probable region of 
parameter space is constrained by the ^H/^He and ^He/^He 
ratio (rather than the best fit to the ^H and '^He fluxes), 
we feel that it is a safe procedure till high precision data 
from PAMELA of AMS-02 are avaflablc. 

The demodulated interstellar (IS) fluxes for ^H and ^He 
are shown in the left panels of Fig. lU whereas the yet 
modulated top-of-atmosphere (TOA) ratios for ^H/^He and 
■^He/^Hc are shown in its right panels. The references for 
the data are given in Table lA.ll 



Appendix B: Cross-sections 

This appendix summarises the production and destruction 
cross-sections employed for the quartet nuclei in this paper. 

B.l. Elastic and inelastic cross-sections 

All reacti on cross-section s are taken from the parametri- 
sations of iTripathi et al.l (jl999t) . but for the pH reaction 
cross-section. The latter is evaluated from cTpp ' = fpp — CTpp, 
where the total and elastic cross-sections are fitteci to the 
data compiled in the PDCj^. Note also that for ''He -|-^ He, 
we had to renormalise iTripathi et aP (jl999t) formulae by a 
factor 0.9 to match the low-energy data. 

Our parametrisations (lines) and the data (symbols) are 
shown i n Fig. [ B.llfor reactio n on H and He. Note that we 
rely on ITripathi et al.l ()1997ll for any other inelastic reac- 
tion. 



|http://pdg. lbl.gov/] 
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Table A.l. References for the quartet data. 



Exp. 



#data 



Year (MV) Ref.' 



Comment 



MASS89 


9 


1989 


1200 


Voyager87 


7 


1987 


360 


Balloon90 


1 


1990 


1200 


Voyager94 


4 


1994 


150 


CAPRICE94 


14 


1994 


600 


IMAX92 


8 


1992 


550 


AMS-01 


10 


1998 


600 


BESS93 


5 


1993 


700 


CAPRICE98 


5 


1998 


700 


MASS89 


5 


1989 


1200 


Voyager87 


7 


1987 


360 


Voyager94 


4 


1994 


150 


IMAX92 


24 


1992 


750 


BESS93 


7 


1993 


700 


BESS94 


5 


1994 


630 


BESS95 


6 


1995 


550 


BESS97 


7 


1997 


491 


BESS98 


7 


1998 


700 


AMS-01 


5 


1998 


600 


IMP3-f4-|-5 


1 


65+67+69 


650 


Balloon73 


1 


1973 


650 


IMP7-fPioneerlO 


3 


72+73 


540 


Balloon77-f Voyager 


14 


1977 


400 


IMP8 


1 


1977 


400 


ISEE3 


1 


78-84 


740 


Balloon84 


1 


1974 


650 


MASS89 


9 


1989 


1200 


Balloon90 


1 


1990 


1200 


IMAX92 


8 


1992 


550 


AMS-01 


4 


1998 


600 


Balloon72 


2 


1972 


400 


IMP7-|-PioneerlO 


2 


72+73 


540 


Balloon73 


2 


1973 


500 


Balloon77+ Voyager 


3 


1977 


400 


Balloon81 


1 


1981 


440 


ISEE3 


2 


78-84 


740 


ISEE3(HIST) 


1 


78 


500 


MASS89 


5 


1989 


1200 


SMILI-I 


12 


1989 


1200 


Voyager87 


1 


1987 


360 


Balloon89 


1 


1989 


1400 


IMAX92 


21 


1992 


750 


SMILI-II 


10 


1991 


1200 


BESS98 


7 


1998 


700 


CAPRICE98 


1 


1998 


700 


AMS-01 


5 


1998 


600 



Webber et al. (1991) 
Seo et al. (1994) 
Bogomolov et al. (_1995) 
Seo fc McDonald (1995) 
Boezio et al. (1999) 
de N olfo ct al. (2000) 
AMS CoU aboration (2002) 
Wang et al. (2002') 
Papini et al. (2004) 

-^He- 
Webber et al. (1991) 
Seo et al. (1994) 
Seo fc McDonald (1995) 
Menn et al. (2000) 
Wang ct al . (2002) 
Myers et al. (2003) 
Myers et al. (2003) 
Myers et al. (2003) 
Myers et al. (2003) 
Aguilar et al. (2011) 

— ^H/^He — 
Hsieh et al. (1971) 
Apparao (1973) 
Teegarden et al. (1975) 
Webber fc Yushak (1983) 
Beatty et al. (1985) 
Kroeger (1986) 
Durgaprasad fc Kunte (1988) 
Webber et al. (1991) 
^ogomolo^_et_aL (1995) 

de Nolfo et al. (2000) 
Aguilar et al. (2011) 

— ^He/^He — 

Webber fc Schofield (1975) 
Teeg arden ct al. ( 1975) 
Leech & Ogallagli^r (1978) 
Webber fc Yushak (1983) 
Jordan (1985) 
Kroeger (1986) 
Mewaldt (1986) 
Webber et al. (1991) 
Beatty et al. (1993) 
Seo et al. (1994) 
jlatano_et_^l^ (1995) 
Menn et al. (2000) 
Ahlen et al. (2000) 
Myers et al. (2003) 
Mocchiutti et al. (2003) 
Aguilar et al. (2011) 



Voyager is at 23 AU 

Voyager is at 56 AU 

Subtraction of H to from Tab 3. 



Voyager is at 23 AU 
Voyager is at 56 AU 



Taken from their Fig. 2 
Taken from their Fig. 2 
Taken from their Fig. 2 
Taken from their Fig. 2 



(discarded in the analysis) 



Re-analysed by Webber et al. (1987 



Re-analysed by Webber et al. (1987) 
(discarded, see ^Vebber^£_al^^987) 



Voyager is at 23 AU 



Supersedes 2ypns,,et_al^ (2003) data 



B.2. Light nuclei production: Nuc + p 

The light nuclei '^He and ■^H are spallativc products of cos- 
mic rays interacting with the interstellar medium (ISM). 
The total secondary flux is obtained from the combina- 
tion of production cross-sections and measured primary 
fluxes. In principle, all nuclei must be considered, but 
the ISM and GCRs are mostly composed of and 
^He, making the reactions involving these species dom- 



inant. For heavier species, their decreasing number is 
balanced by their higher cross-section. In several stud- 



ies (e.g. iRamatv fc Lingenfelterl 119691 : lJung et all Il973bl : 
iBeattvl [T986ir it was found that the CNOcr+ Hism re- 
actions contribute to ^ 30% of the flux above GcV/n 
energies. The reverse reaction Hcr+CNOism mostly pro- 
duces fragments at lower energies, making them irrelevant 
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500 



- Reaction cro?s-sections (on He) 



o ■■ -"He + He 
«- 'He + He 
■- -H + He 



E 



^400 




300 



200 



Fig. B.l. Total inelastic (reaction) cross-section for 
the quartet isotopes. The lines show our parametri- 
sation (see text), the symbols ar e data. Top panel: 
on 



reaction on H with data from 
iHavakawa et all dlOei): lleo et al.l 





Cairns et all (I1964D 
Il967f): iPalevskv et al 



.Griffiths fc HarbisonI (!l96!/): iNichoUs et al . 
Carlson et all (I1973D: S o urkes et al.l (119761): 



lAbleev et al.l (I1977D : iKlem et al.l (ll977D: Tjaros et alJ 




Velichko et al 



lAbdullin et al.l 

Webbeil (Il997f). Bottom p a nel: re action on 
data from iKoepke &: Brownl (|1977t ): iJaros et all (jl978D : 



Blinov et al.l 



19821);^ 

19931 ): iGlagolev e~ 




iTanihata et al.l (|l985[ ). 



for CR studies in the regime > 100 McV/iQ Note that 
is also produced in these reactions, but it decays in ■^Hc 
with a life time (12.2 years) short with respect to the prop- 
agation time. All tritium production is thus assimilated to 
■^He production, but the cross-sections for this fragment are 
provided as well below. 

The energy of the fragments roughly follows a Gaussian 
distribution (e.g. ICucinotta et al.lll993l) . Its impact on the 
secondary f lux w as inspected for the B/C analysis by 
iTsao et all (Il995f ). where an effect < 10% was found, com- 
pared to the straight-ahead approximation, in which the 
kinetic energy per nucleon of the fragment equals that of 
the projectile. The precision sought for the cross-sections 
is driven by the level of precision attained by the CR data 



Solar modulation also ensures that only species created at 
energies > GeV/n matter. 



to analyse. Given the large errors on the existing data, the 
straight-ahead approximation is enough for this analysis. 
However, future high-precision data (e.g. from the AMS-02 
experiment) will probably require a refined description. 

B.2.1. ^Hc + p ^ 3jj^ and ^He 

Recent and illustrative reviews on "^He-f H reaction and the 
produ cti on of light fragm en ts is giv en by B ildstcn et al. 
1990[) : ICucinotta eFall (|1993[) : iBlinov fc Chadeveva 
2008f ). As said earlier, we are only interested in the 
total inclusive production cross-section, not in all the 
possible numerous final states (see, e.g.. Table 3 of 
.Bl inov &: Chadeveva 2008,). We adapt the parametrisation 
of ICucinotta et al-F |i993[ ). which takes into account 
separately the break-up and stripping (for '^He and ^H) 
cross-sections. The former reaction corresponds to the case 
where the helium nucleus breaks up leading to coalescence 
of free nucleons into a new nucleus. The latter happens 
via the pickup reaction where the incident proton tears a 
neutron or a proton off the helium nucleus. Both reaction 
and the total are shown along with the experimental data 
in Fig. |R2l 

The most accurate set of data (upward blue empty 
triangles) ar e from the experiments set up in ITEP and 
LHE JINR (lAladashvill et al.l IT98ll: [Glagolev et al.' TgM 
Abdullin et al .1 11994 summarised in iBlinoy fc Chadeveva 



2008). Their highest energy data point (jClagolev et al 
d993l ) is a conservative estimate as the more or equal 
to 6-prong reactions are not detailed (s ee Table 3 of 
Blinov fc Chadeveval l2008l and Table 4 of iGlagolev et all 



199^ To take into account that possibility, we consider an 
error of a few mb in the plots of Fig. IB. 21 Let us consider 
in turn each product of interest. 



■^Hc production The stripping cross-section data (d and 
• ^He in the fina l stat e) are well fit ted by Eg. (130) of 



ICucinotta et al.l (1 19931). Howev er. the lGriffiths fc HarbisonI 
(|1969D and (jJung et al.lll973a^ are ~ 30% below the other 
data. Actually, for the latter (filled stars) the break-up 
cross-section is above other data, it may be that the end 
products are misreconstructed (in this or the other exper- 
iments). Nevertheless, the sum of the two — which is the 
one that matters — is consistent in all data. Note that we 
slightly modified the break-up cross-section provided by 
ICucinotta et al.l (|l993) to better fit the high-energy data 
points. For the latter, all the data are consistent with 
one another, but for the high precision ITEP data at 200 
MeV/n. 

^H production The stripping cross-section is as for ^He (d 
and "^He in the final state). Th e high-energy b reak-up cross- 
section data (LHE JINR and IWebbedll990b[) are inconsis- 
tent. We have decided to rescale the Webber data, to take 
into accou nt the fact tha. t in his preliminary account of 
the results (|Webbeilll990b[) , the total inelastic cross-section 
i s smaller than that given in a later and updated study 
([Webbeil I1997D . Still, the agreement between the two sets 
is not satisfa ctory. The other high-energy data point is the 
llnnesi (|l957t) experiment, and it suffers large uncertainties 
and maybe systematics (it is for n - j- ^He reaction, and 
the data point is provided by iMeveil (|l972f ) who relied on 
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*He + p ^ -''He 
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Breakup 

Tannenwald 53[Meyer] 
Innes 57 [Meyer] 
Griffiths 69 
Meyer 72 
Nicholls 72 
Jung 73 

izard 77 
Webber 90 

Aladashvili SI, Glagolev 93, Abdullin 94 
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Tannenwald 53 [Meyer] 
Rogers 69 
Griffiths 69 
Meyer 72 
Jung 73 
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E^„ [GeV/n] 
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Innes 57 [Meyer] 
Cairns 64 [Meyer] 
Jung 73 

Webber 90 [rescaled] 



Stripping (d+^e) 
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Meyer 72 
Jung 73 
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Fig. B.2. Inclusive production cross-sections of '^He 
(top), (centre) and (bottom) in "^He-I-H re- 

action. The data ( see t ext fo r details) ar e from 



Tannenwaldl (I1953D : llnnesi (I1957D: ICairris et al I (11964: 
Rogers et all (Il969l) : 'Griffiths & Harbison' (1969"); 'Mcvci) 



(I1972D: lJung et alj (1973a): A la dashv ili ct al.. a981il: 
iWebbeil (Il990bn : iGlagolev et al.l (f!99l : lAbdullin et abl 
(I1994D . 



several assumptions to get it). The ITEP/LHE JINR data 
being the best av ailable, -we have r eplace d the formula for 
the breakup of lCucinotta et all (jl993|) by a form similar 
as that given for ^He, but -where "we changed the parameters 
to fit the high energy points. 
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□ Griffiths 69 




= - X Meyer 72 




: A Blinov 86, Glagolev 93 
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Fig. B.3. other production channels from the less abun- 
dant "^He and the peaked fusion pp reaction (the much 
smaller cross-section is redeemed by a CR flux higher (p 
instead of '^ He). The data are frorn [G riffiths fc Harbison 

mmy . ,gi 



1969 1) : iMevCTl (IT971 : [Blkiov et all 



1991 . 



.Glagolev et al. 



•^H production There is only break-up for the 
ICucinotta et all (jl993l) '^H production. The data are 
in broad agreemeri t -with one another, but for the 
iNicholls et al.l (I1972D point (open p lus). Again, we have 
adapted the ICucinotta et al.l (|l993l) parametrisation to 
better fit the ITEP/LHE JINR data. 

B.2.2. ^Hc + p ^ 2jj (breakup) and p + p ^ ^jj (fusion) 

There are f-wo other channels for producing ^H from light 
nucleus reactions, and they are sho-wn in Fig. IB. 31 along 
-with the data. The first one is from ^He (break-up and 
stripping). The CR flux of the latter is less abundant than 
the "'He flux. With a ratio of ^ 20% at 1 GcV/n (decreasing 
at higher energy) and similar production cross-sections {'^ 
30 — 40 mb), this is expected to contribute by the same 
fraction at GeV/n energies, and then to become negligible 
> 10 GeV/n. The second channel is the ^H coalescence from 
t-wo protons. The cross-section is non-vanishing only for a 
very narro-w energy range. Even if the cross-section is 10 
times smaller than for the other channels, the fact that CR 
protons are ~ 10 times more numerous than ^He makes it 

a significant channel slightly belo-w 1 GeV/n . 

The fitting curves are taken from lMeveJ (|1972D . but "we 
adapted the fit for the ^He-|-p channel to match the t-wo 
high-energy ITEP/LHE JINR data points. 

B.2.3. Proj(^>4) + p ^ 3jj^ grid ^He 

For nuclear fragmentation cross-sections of heavier nuclei, 
the concepts of 'strong' or '-weak' factorisation relies on the 
fact that at high energy enough, the branching of the var- 
ious outgoing particles-production channels becomes inde- 
pendent of the target. This corresponds to the factorisation 
^strong^p^^ ^ ^F^^ cr^<^^i^{p^ F,T) = jPT where 
a{P, F, T) is the fragmentation cross-section for the projec- 
tile P incident upon the targ et T produc i ng the fragment 
F. This is discussed, e.g., in lOlson et all (119831 ). -where it 
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is concluded that although strong factorisation is proba- 
bly violated, weak factorisation seems exact (see also, e.g., 
iMichel et al.l Il995l) . The parametrisation proposed below 
takes advantage of it. 

Several data exist for the production of light isotopes 
from nuclei A> 12 on H (see Tabic [BTT]) . The most com- 
plete sets of data in terms of energy coverage are for the pro- 
jectiles C, N, and O {{A) = 14), the group Mg, Al, and Si 
{{A) = 26), and the group Fe and Ni {{A) = 57). They are 
plotted in Fig. IB. 41 (top panels and bottom left panel). The 
solid lines correspond to an adjustment (by eye), rescaled 

from the ^breakup cross-section (because heavy projectile 
do not give A ~ 3 fragments in the stripping process). As 

'^^He ~ '^^h'^' distinction is made for the fit (^H data 
are scarce and do not influence the conclusions drawn from 
these three groups of nuclei) . The following parametrisation 



aPp-^(i?fe/„, ^p) ^ 7p • fiEk/n,Ap) ■ a-:p;:p"^(i?,/„ 



with 



(B.l) 



fiEk/n,Ap)^< [1.5 GcV/n 



0.8- 



iiEk/n<l-5 GeV/n, 

otherwise; 

(B.2) 

proves to fit well the three groups of data for ener- 
gies greater than a few tens of MeV/n. Thanks to the 
f{E]^/n, ^proj) factor, there is no further energy dependence 
in the 7p factor, so that the latter can be determined from 
the data points at any energy. The bottom right panel of 
Fig. IB. 41 shows the measured mean value and dispcrsiorP^ 
as a function of A. from which we obtain: 



7p 



He 



7p" 



1.3 



Ap 
25 



1.5 



7p 



P = 0.28 Ah"^. 



(B.3) 



The set of formulae (EH]), (|B:2|) . and (|B3)) completely 
define the Proj-|-p production cross-sections for the light 
fragments. 

B.3. Proj(^A>4) +^He-¥^H,^H, and ^He 

Data for T -I- *He where t he target T i s heav ier than p are 
scarce. In a compilation of iDavis et al.l (|1995[ ). the authors 
find that the '^He production scales as A^^^ (based on 4 
data point with At > 7). This is the scaling we employ for 
the ■^H and production as well. 
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Table B.l. References for the light nuclei production cross-sections (references for '^He projectiles are given in Fig. lB.2p . 
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Fig. B.4. Proj(A>4) + P ^ ^H, ^H, and ^He cross-sections for Proj=C,N,0 (top left), Proj=Mg,Al,Si (top right), and 
Proj=Fe,Ni (bottom left). The bottom right correspond to the 7p factor (sec Eg. IB. II and text for explanations). The 
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